Fuente: artículo online

Código

library(sf)
## Linking to GEOS 3.9.3, GDAL 3.5.2, PROJ 8.2.1; sf_use_s2() is TRUE
library(mapview)
## The legacy packages maptools, rgdal, and rgeos, underpinning the sp package,
## which was just loaded, will retire in October 2023.
## Please refer to R-spatial evolution reports for details, especially
## https://r-spatial.org/r/2023/05/15/evolution4.html.
## It may be desirable to make the sf package available;
## package maintainers should consider adding sf to Suggests:.
## The sp package is now running under evolution status 2
##      (status 2 uses the sf package in place of rgdal)
baseurl <- "http://vport.herrera.unt.edu.ar:8894/geoserver/wfs?"
wfs_request <- "request=GetFeature&service=WFS&version=2.0.0&typeName=espacio01:tuc_juris&outputFormat=application/json"

fi_regions_wfs <- paste0(baseurl,wfs_request)
fi_regions <- st_read(fi_regions_wfs)
## Reading layer `OGRGeoJSON' from data source 
##   `http://vport.herrera.unt.edu.ar:8894/geoserver/wfs?request=GetFeature&service=WFS&version=2.0.0&typeName=espacio01:tuc_juris&outputFormat=application/json' 
##   using driver `GeoJSON'
## Simple feature collection with 121 features and 4 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -66.17672 ymin: -28.02735 xmax: -64.48182 ymax: -26.05595
## Geodetic CRS:  WGS 84
mapview(fi_regions, label = fi_regions$name, color = "darkgreen", col.regions = "green", alpha.regions = .05)

Explicación

Con gusto te explico el código!

Cargando las librerías necesarias

library(sf)
library(mapview)

library(sf): Esta línea carga la librería sf, que es fundamental para trabajar con datos espaciales en R. sf nos permite leer, manipular y analizar datos geográficos de manera eficiente.

library(mapview): Cargamos la librería mapview para crear mapas interactivos de manera sencilla. Con esta librería, podremos visualizar nuestros datos espaciales en un mapa web. Construyendo la URL del servicio WFS

Importación del WFS (Localización+Ubicación+unión de la solicitud)

baseurl <- "http://vport.herrera.unt.edu.ar:8894/geoserver/wfs?"
wfs_request <- "request=GetFeature&service=WFS&version=2.0.0&typeName=espacio01:tuc_juris&outputFormat=application/json"

fi_regions_wfs <- paste0(baseurl,wfs_request)

baseurl: Aquí se define la dirección base del servidor WFS. Un servidor WFS (Web Feature Service) es un servicio web que permite acceder a datos geográficos en formato vectorial.

wfs_request: Esta parte especifica la solicitud que se le hará al servidor WFS. En este caso, estamos solicitando un conjunto de características separados por & request=GetFeature: Indica que se solicita obtener datos de una capa (feature) del servicio WFS. service=WFS: Especifica que se trata de un servicio WFS. version=2.0.0: Define la versión del servicio WFS que se está utilizando (en este caso, la versión 2.0.0). typeName=espacio01:tuc_juris: Define el nombre completo de la capa que se desea descargar. (“espacio01:tuc_juris”) outputFormat=application/json: Indica que se desea recibir los datos en formato JSON.

fi_regions_wfs: Finalmente, se combina la dirección base con la solicitud para formar la URL completa del servicio WFS.

Leyendo los datos del servicio WFS (crea un archivo sf, a partir de la solicitud)

fi_regions <- st_read(fi_regions_wfs)
## Reading layer `OGRGeoJSON' from data source 
##   `http://vport.herrera.unt.edu.ar:8894/geoserver/wfs?request=GetFeature&service=WFS&version=2.0.0&typeName=espacio01:tuc_juris&outputFormat=application/json' 
##   using driver `GeoJSON'
## Simple feature collection with 121 features and 4 fields
## Geometry type: MULTIPOLYGON
## Dimension:     XY
## Bounding box:  xmin: -66.17672 ymin: -28.02735 xmax: -64.48182 ymax: -26.05595
## Geodetic CRS:  WGS 84

st_read(fi_regions_wfs): Esta función de la librería sf se utiliza para leer los datos espaciales desde la URL del servicio WFS. Los datos se leen y se almacenan en un objeto de tipo sf, que es la forma estándar en R para representar datos espaciales.

Visualizando los datos en un mapa

mapview(fi_regions, label = fi_regions$name, color = "darkgreen", col.regions = "green", alpha.regions = .05)

mapview(fi_regions): Esta función crea un mapa interactivo utilizando los datos espaciales almacenados en fi_regions. label = fi_regions$name: Cada polígono en el mapa se etiquetará con el valor correspondiente de la columna “name” en los datos. color = “darkgreen”, col.regions = “green”, alpha.regions = .05: Estos argumentos se utilizan para personalizar la apariencia del mapa, definiendo el color del borde de los polígonos, el color de relleno y la transparencia del relleno, respectivamente.

En resumen

Este código completo realiza las siguientes tareas:

Conecta a un servidor WFS: Se establece una conexión a un servicio web que proporciona datos geográficos en formato vectorial.

Descarga datos: Se descargan los datos de una capa específica (en este caso, “espacio01:tuc_juris”) del servidor WFS.

Crea un objeto espacial: Los datos descargados se convierten en un objeto espacial en R, lo que permite realizar análisis espaciales.

Visualiza los datos: Los datos espaciales se muestran en un mapa interactivo, donde cada polígono está etiquetado con un nombre y tiene un estilo de color específico.

En pocas palabras: este código te permite visualizar datos geográficos de un servidor remoto en un mapa interactivo utilizando R.